Band structure analysis of the conduction-band mass anisotropy in 6H and 4H SiC 
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The band structures of 6H SiC and 4H SiC calculated by 
means of the full-potential linear-mufhn-tin orbital method 
are used to determine the effective mass tensors for their 
conduction-band minima. The results are shown to be consis- 
tent with recent optically detected cyclotron resonance mea- 
surements which find the ratio of cyclotron masses for B 1 c 
to B || c to be larger (smaller) than unity for the 6H (4H) 
polytype. However, contrary to previous suggestions, appre- 
ciable anisotropies in the c-plane are found. For 6H SiC, a 
strong dependencies on band-filling is predicted because of 
the occurence of a double well minimum along the ML-axis. 
The calculated mass tensors for 3C and 2H are also reported. 

PACS numbers: 71.25.Jd 71.25.Tn 



The technological potential of SiC for high-tempera- 
ture, high-power and high-frequency electronic devices as 
well as the intrinsic scientific interest of the polytypism 
of SiC have been recognized for several decades. Only 
recently, however, has it become possible to obtain elec- 
tronic grade quality single crystal single polytype mate- 
rial. Recently, it was found that 4H SiC exhibits a much 
smaller (and opposite) anisotropy of the electrarumobil- 
ity with respect to the c-axis than does 6H SiCmj It has 
been unclear up to now whether this could be attributed 
to a corresponding anisotropy in the effective mass ten- 
sors. Rather conflicting results have been reported on the 
effective masses as determined by a fitting of LB|-|absorp- 
tion spectra of the excited shallow donor statesa'Qand by 
optically detected cyclotron resonance (ODCR).LrQ 

In this letter, calculated band structures of 6H SiC and 
4H SiC are used to determine the effective mass tensors at 
the conduction-band minimum. Our results are in good 
agreement with the ODCR results on the anisotropies 
with respect to the c-axis for 6H and 4H SiC masses and 
hence support a mass-anisotropy as the cause for the mo- 
bility anisotropies. Our calculations also predict appre- 
ciable anisotropies of the mass tensors in the c-plane. In 
addition, they predict a strong dependency of the masses 
on band-filling in the case of 6H SiC due to a double- 
well-like conduction-band minimum. 

Our band-structure calculations were performed 
within the density functional theory in the local den- 
sity approximation (LDA)El usingj-the full-potential (FP) 
linear muffin-tin orbital method.Bt] They differ in some 
small, but for the .present purpose relevant details, from 
our earlier resultal3_which used the atomic sphere ap- 
proximation (ASA) O The previous calculations were al- 
ready found to provide a satisfactory account of the UV- 



reflectivity spectra.0 In the present work, we performed 
calculations on a fine mesh of k-points in order to exam- 
ine the details of the energy-band surfaces close to the 
conduction-band minima. 

In 6H SiC, the minimum is found about halfway be- 
tween M andpX along the ML = t/-axis of the 6H- 
Brillouin zone.tJ This leads to a 6- valley model in which 
the valleys occur in 3 pairs with each pair being separated 
by a small barrier (of 9.5 meV) at the three equivalent 
M points. The lowest band along this line can be well 
approximated by a fourth order polynomial in \k z \, with 
z chosen along c. In the directions parallel to Mr (or 
[1010]) and MK (or [2110]), respectively called the x- 
and y-directions, the bands are found to behave parabol- 
ically near the minimum with the masses there almost 
equal to those at M . The cyclotron masses for orbits not 
in the c-plane are calculated from m = (h 2 /2n)dA/dE, 
with A the cross-sectional area of the extremal orbit on 
the constant energy surface calculated numerically from 
the fitted energy bands. By varying the energy E of the 
chosen constant energy surface, we obtain the variation 
of the mass with conduction bandjiJling. 

In 4H SiC the minimum is foundEJ at M and the bands 
are parabolic in the three symmetry directions of the 
crystal up to at least 300 meV above the minimum. In 
fact, there are two bands at M within less than 100 meV 
of each other whose energetic position was found to be 
sensitive to the details of the calculation. Specifically 
these bands were interchanged in the ASA calculation. 
Apart from this interchange, the results from FP and 
ASA calculations were in very good agreement (includ- 
ing the masses of each of these bands). 

The first three rows in Table [| give the principal com- 
ponents of the mass tensor at the ML minimum in 6H 
SiC and at the M minima in 4H and 2H SiC. Clearly, 
these are fully anisotropic (i.e. mi ^ rri2 ^ 7713 with 
1,2,3 corresponding to x,y,z.) This is consistent with 
the Ci v symmetry for points on the ML axis. 2H SiC 
is included although its lowest conduction band mini- 
mum is at K (in the equatorial plane at the corner of 
the hexagon.) It can be seen that the 2H masses at M 
are more similar to those of 6H than those of 4H. Similar 
closer similarity between 6H and 2H than between 4H 
and 2H was previously also found in other aspects of the 
band structure.lia The mass tensor at K is rotationally 
invariant about the c-axis as required by the C$ v sym- 
metry of the AT-point and has the values: m || c = 0.27 
and m _L c = 0.45. 

The cyclotron effective mass associated with a single 
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elliposid is: 

/ m 1 m 2 m 3 
y niibf + m 2 02 + TO3O3 

where hi are the direction cosines of B with the principal 
axes of the mass tensor. To compare our results with 
the ODCR measurements it is necessary to consider the 
combined effect of the different valleys which for a general 
orientation of the magnetic field B become inequivalent. 
For B along the c-axis, the 3 valleys for 4H (or double val- 
leys for 6H) at the different M-points are equivalently ori- 
ented; consequently only one peak is observed in ODCR 
with an effective mass m* ± = ^Jm^m^. Our values for this 
effective mass are given in the sixth row of Table [j| and are 
in good agreement with the ODCR values of 0.42 ± 0.02 
in both 4H and 6H SiC. The slight underestimate of the 
calculated masses for 4H is partially due to the absence of 
phonon renormalization contributions and may also par- 
tially be due to the LDA. In this connection we note a 
similar situation in 3C-SiC. There, our calculated values 
for the electron mass are mi = 0.63, mt ^ 0.23, while the 
values measured by cyclotron resonancet£l are mi = 0.677 
and mt = 0.247. (Here longitudinal (I) and transverse (t) 
refer to the directions parallel and perpendicular to the 
cubic r — X axis for the A-minimum.) 

For B in the c-plane, the ellipsoids at the various M 
or ML points are generally oriented differently with re- 
spect to the field and should give rise to three peaks in 
ODCR. For the [1010] and [2110] directions, two of the 
peaks coincide. The variation of the masses in the c- 
plane is shown in Fig. |as a function of azimuthal angle 
<j> measured from the [1010] direction. 

The published ODCR measurements using microwave 
frequencies of 9.235 GHz in the X-band do not report 
anisotropy in the c-plane and were analysed as if the 
valleys were ellipsoids of revolution with a m± and TO||, 
where _L and || refer to the c-axis. Fig. || shows simu- 
lated ODCR microwave absorption spectra based on the 
calculated masses for 4H SiC and 6H SiC for a number 
of important directions. The lot values (microwave fre- 
quency times lifetime) used are cloaa,to those deduced 
from the experiments of Son et ai£m. Because of the 
width of the ODCR signals from the differently oriented 
ellipsoids, the in-plane sum peak is seen to show little 
variation in position. This explains why the ODCR spec- 
tra appear to give an isotropic in-plane mass in spite of 
the fact that each ellipsoid is strongly anisotropic in the 
plane. When the effective in-plane masses are equated to 
^m^m|, as is appropriate for a model with in-plane ro- 

tationally invariant ellipsoids, the effective m|'s in Table 
I are obtained. 

Son et aln also reported an ODCR signal using a 
higher microwave frequency in the Q-band (y — 35 GHz). 
With the correspondingly higher lot = 9.1 value, our cal- 
culations predict that the mass anisotropy in the c-plane 
should be clearly resolved for 4H, as can be seen in Fig. [|. 
Very recently, such in-plane anisotropy was confirmed by 



Hofmann et aiJ13. The mass values deduced from these 
new experiments (mi = 0.58, m 2 = 0.31 and 7713 = 0.33) 
are in good agreement with our calculated values. 

For 6H SiC, the strong non-parabolicity resulting from 
the double well character of the minima along the ML 
axis leads to a variation of the masses with energy above 
the minimum (E — E c ) (or, band filling). Fig. || shows 
the variation with E — E c of v /mini3, \Jm2mj, and the 
effective mjj deduced from the calculated broad peak in 
microwave absorption for B _L c in the manner explained 
above. The calculated m^ is seen to gradually increase 
with band filling up to the barrier, at which point there 
is a discontinuity because the orbits then encircle the 
double minimum centered at M instead of the separate 
valleys. At higher energies mj| decreases with increas- 
ing energy. In the degenerate limit, a Fermi energy just 
at the barrier corresponds to a carrier concentration of 
the order of 10 19 cm -3 . While carrier concentrations of 
this magnitude are certainly possible n the carrier concen- 
trations in the ODCR measurementsQ (at 6K) which are 
produced by optical pumping are expected to be much 
lower. Based an. the laser intensity employed and the 
carrier lifetimesEj we estimate concentrations of the or- 
der of 10 15 cm~ 3 . On the other hand, the high microwave 
powers used (> 20 mW) may lead to a non-equilibrium 
hot carrier distribution and hence the experiment may 
conceivably probe masses higher in the band. 

As can be seen in Table Q, the two calculated 4H masses 
and m*j_ for 6H SiC arc in excellent agreement with ex- 
periment. Also consistent with experiment, mjj in 6H 
SiC is relatively large. Nevertheless, this mass is 40 % 
smaller than experiment. Several factors, including the 
non-parabolicity noted above could bear on this discrep- 
ancy. The polaron renormalization mass increase could 
be substantial for this large mass because the polaron 
coupling constant increases with mass. Electron-electron 
interaction corrections to the LDA may also alter the 
masses, but are unlikely to be significant in view of the 
good agreement for the other masses. Finally, we note 
that the experimental uncertainty on this mass is fairly 
large because it depends on fitting a broad resonance. 
Considering the above factors, the agreement between 
the calculated and measured values for 6H mj| may be 
called fair. 

In conclusion, the agreement between our calculated 
effective mass tensors and ODCR measurements is ex- 
cellent for 4H SiC and semiquantitative for 6H SiC. Al- 
though we found the effective mass tensors to be fully 
anisotropic in both polytypes, we explain why the recent 
X-band ODCR measurements did not resolve the super- 
position of the fairly broad resonant peaks for B 1 c. 
That lead to the erroneous suggestion that the in-plane 
mass is isotropic. In agreement with the ODCR stud- 
ies, we find the ratio m|/m^ < 1 for 4H and > 1 for 
6H SiC. The in-plane anisotropy is predicted by our cal- 
culations to be observable with higher lot values as was 
confirmed very recently for 4H SiC by ODCR measure- 
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merits using a higher microwave frequency.113 We fur- 
ther predict a considerable dependence of mj in 6H SiC 
on band filling due to the double well character of the 
lowest conduction band. This could influence transport 
properties and may also be detectable in future ODCR 
measurements with higher band fillings. The manner in 
which the effective masses enter the mobility anisotropy 
is far from trivial because the different scattering pro- 
cesses that could come in have different functional de- 
pendence on the masses. Nevertheless, whatever scatter- 
ing mechanisms are involved, the mobility is expected to 
decrease when the mass increases. The calculated c-axis 
anisotropics of the mass tensor thus appear to be the 
basis of the observed large Hall mobility anisotropics in 
6H and small and reversed one in 4H. However, further 
detailed work involving the effects of different scatter- 
ing mechanisms will be required to quantitatively explain 
their magnitude and temperature dependence. Finally, it 
is clear that these rather unusual conduction-band min- 
ima call for an extension of effective mass theory before 
one can confidently deduce effective masses from donor 
excitation IR spectra. 
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TABLE I. Principal components of the conduction-band 
effective mass tensor (in units of the free electron mass) for 
6H (ML-minimum), 4H and 2H SiC (at M- minimum) and 
derived quantities. 
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FIG. 1. Cyclotron masses in c-plane for 6H and 4H SiC as 
function of azimuthal angle cj>. <j) — corresponds to [1010] 
and <j> = 30° to [2ll0]. The line type corresponds to the 
individual contributions to the ODCR spectra shown in Fig. 
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FIG. 2. Simulated X-band (9.235 GHz) ODCR spectra for 
4H SiC and 6H SiC. Solid line, B || [0001]; dashed line, 
B || [1010]; dotted line, B || [2110] 
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FIG. 4. Calculated effective masses in 6H SiC as a function 
of conduction band filling. The six-valley region below the 
barrier is shown separately from the three- valley region above. 
m xz and m yz are defined in Table I. mj] is a pseudo-effective 
mass parallel to c appearing in the isotropic in-plane model 
discussed in the text. 



FIG. 3. Simulated Q-band (35 GHz) ODCR spectra for 4H 
SiC using lot = 9.1 for various orientations of the magnetic 
field. The dashed, dotted and dot-dashed lines correspond to 
the individual effective mass contributions shown in Fig. ^ 
and the solid lines to their sums. For clarity, slight off-sets 
are included where two identical peaks overlap. 
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